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ABSTRACT 

We study the radio spectrum of PSR B 1259-63 orbiting around the Be star LS 2883 and 
show that the shape of the spectrum depends on the orbital phase. At frequencies below 3 GHz 
PSR B1259-63 flux densities are lower when measured near the periastron passage than those 
measured far from periastron. We suggest that an interaction of the radio waves with the Be 
star environment accounts for this effect. While it is quite natural to explain the pulsar eclipse 
by the presence of an equatorial disk around LS 2883, this disk alone cannot be responsible 
for the observed spectral evolution of PSR B1259-63 and we, therefore, propose a qualitative 
model which explains this evolution. We consider two mechanisms that might influence the 
observed radio emission: free-free absorption and cyclotron resonance. We believe that this 
binary system can hold the clue to the understanding of gigahertz -peaked spectra of pulsars. 

Key words: (stars:) pulsars: general - radiation mechanism: non-thermal - pulsars: individ- 
ual: B1259-63, B1054-62, J1809-1917,B1822-22, B1823-13, B1828-11, - stars: winds, 
outflows 



1 INTRODUCTION 

Generally, the observed radio spectra of most pulsars can be mod- 
elled as a power la w with negative spectral indices of about -1.8 
dMaron et alj|2000b . If a pulsar can be observed at frequencies low 
enough (i.e. 100 - 600 MHz) it may also show a low freq uency 
turn-over in its spectru m dSieber]|l973l ; lMalofeev et al.ll 1994b . 

On the other hand lLorimer et al J d 1995b mentioned three pul- 
sars which have positive spectral indices in the frequency range 
300 - 1600 MHz. Later iMaron et alj d200fjb reexamined spectra 
of these pulsars taking into account the data obtained at higher 
frequencies (above 1.6 GHz) and consequently were the first to 
demonstrate a possible existence of spectra with turn-over at high 
freq uencies, about 1 GHz. Motivated by this unusual spectral fea- 
ture IffijajT&^aroiil |2004) selected several pulsars w hich showed 
a decr easing flux density at frequencies below 1 GHz. Kija k et all 
d2007l) used multifrequency flux measurements for these candidate 
pulsars and presented the first direct evidence of a high frequency 
turn-over. A frequency at which such a spectrum shows the maxi- 
mum flux was called the peak frequency. 

Based on their observations of these pulsars lKiiak et al.1d201 lb 
provided a definite evidence for a new type of pulsar radio spec- 
tra. These spectra show the maximum flux above 1 GHz, while at 
higher frequencies the spectra look like a typical pulsar spectrum. 
At lower frequencies (below 1 GHz) the observed flux decreases, 
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showing a positive spectral index. They called these objects the 
gigahertz-peaked spectra (GPS) pulsars. 

iKiiak et all d201 lb also indicated that the GPS pulsars are rel- 
atively young objects and they usually adjoin such interesting en- 
vironments as HII regions or compact pulsar wind nebulae (PWN). 
Additionally they seem to be coincident with the known but some- 
times unidentified X-ray sources from 3rd EGRET Catalogue or 
HESS observations. We can assume that the GPS appearance owes 
to the environmental conditions around the neutron stars rather than 
to the radio emission mechanis m. The issue of our interest PSR 
B 1259-63 was also inscribed bv lLorimer eTail dl995h in the list of 
pulsars with positive spectral indices. Therefore, B 1259-63 seems 
a natural candidate to be classified as the GPS pulsar. This pulsar 
was discovered by I Johnston et alj d 1992b in 1990 during a large- 
scale high frequency (1.5 GHz) survey of the Galactic plane. It is 
important to mention that PSR B 1259-63 is the only known radio 
pulsar orbiting a massive main-sequence Be star. 

PSR B 1259-63 has a short period of 48 ms and a characteris- 
tic age of 330 kyr. Its average DM is about 147 pc cirr 3 and the cor- 
responding distance is about 2.75 kpc. The companion star LS 2883 
is a 10-mag massive Be star with a mass of about 10 Mq and a ra- 
dius of 6 Rq. Let us note that Be stars are generally believed to 
have a hot tenuous polar w ind and a cooler high den sity equatorial 
disk. I Johnston et alj < Tl994l) and lMelatos et al.ldl995h suggested the 
presence of the disk around LS 2883. Its density is quite high near 
the star and falls off as a power law with the distance from the star. 
The disk is likely to be highly tilted with respect to the pulsar or- 
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Table 1. The results of our spectral fits to the known GPS pulsars. The 



GPS pulsars 



values of peak frequencies are taken from Kijak et alJfeOlll 
timated using different spectral fit methods (see text for details). In the 
first column (PSR) numbers correspond to the following pulsars B 1054-62, 
J1809-1917, B1822-14, B1823-13, B1828-11 consequently. 



PSR 


a 


b 


c 


x 2 


v K 
p 


Vp 


1 


-7.28±1.98 


-1.95±0.50 


1.94±0.07 


4.40 


1.0 


0.73 


2 


-3.96+1.14 


2.80±0.84 


-0.10±0.09 


2.91 


1.7 


2.25 


3 


-1.30±0.29 


0.40±0.20 


0.41±0.03 


1.59 


1.4 


1.43 


4 


-1.18±0.22 


0.52±0.17 


0.56±0.04 


2.35 


1.6 


1.66 


5 


-2.09±0.24 


-0.55±0.11 


0.12±0.03 


1.69 


1.2 


0.74 



bital plane, and PSR B 1259-63 is eclipsed for about 35 days as it 
goes behind the disk. The pulsar has a long orbital period of 1237 
days and a large eccentricity of 0.87 with a projected semi-major 
axis, a sin ( of 2.6 AU. 

The PSR B 1259-63/LS 2883 is a unique binary system, which 
emits unpulsed non-thermal emission over a wide range of frequen- 
cies ranging between radio and y-rays and its flux varies with or- 
bital phase. Chandra observed this system in 2009 (shortly after 
the apastron passage) and faile d to detect any pul sation with a pe- 
riod equal to that of the pulsar dPavlov et l|). The broadband 
emission from this system is generally believed to be produced by 
the ultra-relativistic electrons accelerated in the bow-shock that ap- 
pears while the relativistic pulsar wind (PW) collides with the dense 
non-relativistic outflo w of the companio n star, producing the Pulsar 
Wind Nebula (PWN) dKong et alj2oTTh . 

In this paper we analyse the available published data of PSR 
B 1259-63 radio flux measured during three periastron passages (in 
1997, 2001 and 2004) and present detailed study of the pulsar spec- 
trum. 



2 PULSARS WITH THE GPS 

To compare various shapes of the B 1259-63 spectrum observed 
at different orbital phases with that of the GPS pulsars we use 
the same fitting method. T his method was previously used by 
iKuzmin & Losovskvl d200ll) for pulsar spectra with turn-over at 
low frequencies. Using the data of flux measurements with er- 
rors for several GPS pulsars (or at least showing a s ignatu re of 
this unusual spec tral feature) taken from iKiiak et al.l d2007b and 
iKiiaketalJfeoilh we analysed the spectra using a function 



F(v) = 10 



here x = log 1( 



(i) 



We carried out fitting of the data using an implementation of the 
nonlinear least-squares Marquardt-Levenberg algorithm. 

In Table 1 we present the values of fitting parameters a, b and 
c, the reduced x 2 , the peak frequency calculated using th e method 
described above and vf, are taken from lKiiak et"al] feOllh . The fit- 
ted spectra are presented in Fig.Q] 



3 BE STAR/PULSAR BINARY AND PSR B1259-63 
SPECTRUM EVOLUTION 

To demonstrate the unusual spectral behaviour of B 1259-63 let us 
fit its spectra applying the typical power law method. Using the 
available measurements of the pulsed flux dManchester & Johnstonl 
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Figure 1. The spectra of five known G PS pulsars. The flux measurements 
are taken from Kijak et al. (2007, 201 1), and the curves represent our fits to 
the data points (see Table 1). 



PSRB1259-63, J1302-6350 




Frequency (GHz) 

Figure 2. The average spectra of PSR B 1259-63 obtained for three differ- 
ent periastron passages (1997, 2000 and 2 004). Three black dots represen t 
the flux measurements obtained in 1993 by Manchester & Johnston (1995). 



1995 1: IJohnston et ail 1 19991 : IConnors et all 120021 : Ijohnston et all 
2005), we have constructed spectra of PSR B 1259-63 for few dif- 
ferent orbital epochs (see Fig. [2jl. The obtained values of spectral 
indices differ signi ficantly from the typical value of about -1.8 
dMaron et alJl200Cj) . However, Fig. [2] shows that this approxima- 
tion is apparently inaccurate as the spectra clearly show the peaked 
frequency for the corresponding periastron passages (see Fig.Q]for 
comparison). Therefore, we have decided to study variation of the 
observed radio flux with orbital phases. 

Using the same database (1997, 2000 and 2004) we have cal- 
culated an average flux density for four given frequencies and for 
each chosen interval of orbital phases. Fig.[3]shows the spectra for 
various orbital phase ranges, determined by the chosen days prior 
to or past the periastron passage. It is clear that the flux at the given 
frequency apparently changes with orbital phases. When the pul- 
sar is close to periastron, the flux generally decreases at all ob- 
served frequencies, but its drastic decrease is observed at the low- 
est frequency. The spectrum for the furthest orbital epoch during 
of which the pulsar was observed so far (i.e. from 113 to 186 days 
past periastron; see the green thickened segment in Fig. [4] bottom 
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Table 2. The results of our spectral fits to the PSR B 1259-63 flux densities, 
averaged over the given orbital epochs, determined by the chosen days prior 
to or past the periastron passage (see also Figs 3 and 4). There are no errors 
in the fit for the orbital phase corresponding to day -18 , because in this 
degenerate case, all errors are zero by definition (numbers of data points 
and parameters are the same, see Fig. [5] panel a). 
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days 
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x 1 


v p 












(GHz) 


-60 : -40 


-1.70±0.46 


1.38±0.54 


0.27±0.16 


0.74 


2.54 


-24 : -21 


-2.10±0.22 


1.76±0.23 


0.10±0.05 


0.10 


2.62 


-18 


-2.77 


2.84 


-0.89 




3.26 


16 : 20 


-4.59±0.03 


4.84±0.02 


-0.71±0.01 


0.01 


3.37 


21 : 24 


-1.72±0.05 


1.51±0.06 


0.31±0.01 


0.01 


2.75 


27 : 55 


-0.88±0.35 


0.34±0.38 


0.55±0.09 


1.49 


1.57 


63 : 94 


-1.48±0.10 


0.55+0.11 


0.56±0.03 


0.04 


1.53 



panel) is consistent with a typical pulsar spectrum described by a 
power-law function (see Fig.[3]panel c and Fig.|4]left panel). On the 
other hand, the spectral index for this orbital epoch is of about -0.8 
which is less than the average value (—1.8) for the pulsar population 
and one can argue that it rather resembles a broken type spectrum 
dMaron et al]|2000h . Thus, it cannot be ruled out that this spectrum 
is still affected by the stellar wind. Observations of the pulsar near 
the apastron passage can clarify this problem. 

To analyse the spectrum evolution in more details we have 
used the method of fitting described in Section 2. Consequently, we 
have calculated fitting parameters for function (1) and constructed 
radio spectra for the chosen orbital epochs. The fitting parameters, 
the reduced^ 2 and the calculated peak frequencies are presented in 
Table 2, while the shapes of the spectrum are presented in Fig. [4] 
left and right panels. Let us note that we have used an additional 
black dashed line in the right panel of Fig.|4]to indicate a fit with an 
additional predicted data point which corresponds to the probable 
flux at frequency 2.6 GHz as there are no observations for this par- 
ticular orbital phase (+18 day) at this frequency. The solid black 
line corresponds to the fit obtained using only three data points 
from Fig. [3] On the other hand one can see that the flux at 2.6 
GHz does not change significantly for all orbital phases (see Fig. 
[3). Thus, in order to use the same kind of fitting for all chosen or- 
bital epochs we have calculated the extrapolated flux at 2.6 GHz 
based on a power law fitting of two data points at higher frequen- 
cies (Fig. [3] panel a). 

The small (much less than unity) values of the reduced x 2 m 
Table 2 are caused by rather large uncertainties in the flux mea- 
surements for four orbital epochs (see Fig.3 for comparison). Let 
us note, that x 2 is especially small (of about 0.01) for +16 : 20 
and +21 : 24 epochs when the relative uncertainties of the mea- 
sured flux are quite large (see Fig. 3b). On the other hand, in the 
case of epochs when the uncertainties (for all frequencies) are rel- 
atively small (see the blue filled dots in Fig. 3a and the blue empty 
diamonds in Fig. 3c) the value of x 1 is close to unity. We hope to 
reduce the uncertainties significantly by gathering more flux mea- 
surements during these orbital epochs. 

To clarify the spectrum evolution we have presented the corre- 
sponding orbital epochs in the bottom panel of Fig. [4] Our analysis 
shows that the shape of the PSR B 1259-63 spectrum depends on 
the orbital phase and therefore it definitely undergoes evolution. It 
should be underlined that the peak frequency also depends on the 
pulsar orbital phase. Comparing Figs Q] and [4] (left and right pan- 
els) we can conclude that the spectra of B 1259-63 resemble those 
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Figure 3. The spectra of PSR B 1259-63 averaged over each orbital epoch, 
from 60 days prior to periastron (panel a) up to 186 days after it (panel c). 



of the GPS pulsars. Moreover, we can see that the shapes of the 
B 1259-63 spectrum at different orbital phase intervals mimic those 
of various gigahertz-peaked spectra. 



4 DISCUSSION 

Multi-wavelength observations showed that some pulsars with 
turn-over in the spectra at high frequencies have very interest- 
ing environments. For exam ple, PSR B 1054-62 l ies behind or 
within a dense HII region dKoribalski etal] Il995h . while PSR 
B 1823- 1 3 surroundings appea r to show some peculiar properties 
in radi o (Gaensl eret al.|[2003T) . as well as in X-ray observations 
jPavlov et al.ll2008l) . which may indicate the existence of a com- 
pact Pulsar Wind Nebula (PWN). The same holds for the GPS pul- 
sar J1809-1917 I Kargaltse v& Pavlov! 1200^ . This could suggest 
that the phenomenon of the turn-over at high frequencies is associ- 
ated with environmental conditions around the neutron stars rather 
than with the radio emission mechanism. Because of more or less 
stable environmental conditions around the GPS pulsars the shapes 
of their spectra do not vary on observable timescales. But in the 
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case of B 1259-63, the pulsar environment considerably changes 
due to high orbital eccentricity while it goes through various orbital 
phases. The distance of the pulsar from the Be star varies from 4.8 
AU in apastron to 0.34 AU in periastron (by a factor of 14). Be stars 
have a strong stellar wind and possibly a strong magnetic field and 
both the disk density and the magnetic field decrease as the dis- 
tance from the star increases. Thus, it is natural to expect that the 
inhomogeneous environment alters the spectrum of B 1259-63 in 
different ways. 

At present it is difficult to construct a detailed theory of the 
spectra evolution as we lack observational data especially near 
apastron (see Fig.[4]l. These data must be important because the pul- 
sar spectrum is the least affected by the environment near apastron 
and this spectrum could be used as a reference spectrum. However, 
from the available limited observations we can still deduce general 
conclusions about the main factors that influence variations of the 
observed spectra. 

Generally, both the hot stellar wind (by means of ther- 
mal absorption) and t he magnetic field may provide absorp- 
tion a t low frequencies (Sieber 1973; Khechinas hvili & Me likidze 
1 19971 : iKhecMnashyiiT et al. 200fJ. While the eclipse itself can 
be naturally explain ed by free-free absorption in the stellar disk 
(Mel atos etafll 19951) . the disk alone is not enough to explain the 
spectra evolution. 

The sign of the GPS feature first appears as early as about 
60 days prior to periastron and it is still observed during at least 
90 days after periastron. It is clear that the stellar disk cannot ex- 
tend that far. It seems that the polar wind is the only factor that 
can affect the spectrum of the pulsar. Preliminary estimations show 
that free-free absorption due to the stellar wind can be responsi- 
ble for this effect before the eclipse as well as after it (see Fig. [4j. 
The optical depth Tff of fr ee-free absorption can be expressed as 
dRvbicki & LightmarJI 19791) : 



T ff = 0.4 x t; 3/2 v GHi 



n.dh 



(2) 



where n e is the electron density in cm 4 , is the temperature in 
10 3 K, v G Hz is the frequency in GHz and / AU is the distance in AU. 
The exact value of T// depends on the geometry of the binary sys- 
tem as well as on the density and temperature of the stellar wind. 
One can note that the shapes of spectra are not fully symmetric with 
respect to the periastron point. To some extent, the spectrum plot- 
ted by the green curve in the left panel of Fig. [4] can be used as a 
reference spectrum. This spectrum is obtained while the pulsar is in 
the orbital epoch designated by the thickened green segment in the 
bottom panel of Fig. [4] It is the furthest orbital phase at which the 
pulsar is observed so far, but the spectrum still does not look exactly 
like a typical pulsar spectrum. Indeed, the waves emitted before the 
eclipse travel longer distance through the stellar wind, thus they are 
stronger attenuated by free-free absorption, than the waves emitted 
after the eclipse (see Fig.[4]l . This is in agreem ent with DM and RM 
variations presented by Mel atos et a 1. (1995). The most significant 
change of the spectrum shape occurs during a couple of days just 
before the eclipse as well as immediately afterwards (the thickened 
black segments in the bottom panel, Fig. |4). One can see that the 
difference between the red and black curves on the left as well as on 
the right panel of Fig.|4]is much more noticeable than the difference 
between the red and blue curves in the same panels. 

In addition, we have shown (see Table 2) that the peak fre- 
quency (as it follows from the fitting procedure) also depends 
on the orbital phase. This effect suggests that the peak frequency 
varies with the changes of the pulsar environment. We argue that 



such behaviour can be explained by the radio-wave absorption in 
the magnetic field associated with the disk. In this case the ra- 
dio waves should pass through the 'magnetosphere' of the disk. 
Magnetic field lines just above the stellar disk are populated by 
the electrons and positrons of the pulsar wind. Such kind of en- 
vironment affects the pulsar spectra in the same way as it does 
in the case of those eclipsing puls ars whose eclipse duration 
depends on the observed frequen cy dKhechinashvili & Melikidze 
ll997l : lKhechinashvili et aT]|2000h . The frequency of waves which 
are affected by the cyclotron dumping can be expressed in the fol- 
lowing way: 



v GHz * 2.8 x 10~ 3 



S 



y p (l - cos0)' 



(3) 



where Bo is the value of the magnetic field associated with the disk, 
y p is the Lorentz factor of the secondary electron-positron pairs of 
the pulsar wind and 6 is the angle between the wave vector and 
magnetic field direction. 

Finally let us underline that the PSR B1259-63/LS 2883 bi- 
nary system can hold the clue to the understanding of the GPS pul- 
sars. As we mentioned in Section 3 the spectrum of B 1259-63 at 
the various orbital phases mimics that of the pulsars with GPS. 
Thus, we can conclude that the GPS feature should be caused 
by some external factors rather than by the emission mechanism. 
On the other hand, the GPS pulsars are isolated radio pulsars 
and therefore, we cannot draw a direct analogy between the PSR 
B1259-63/LS 2883 system and the the GPS pulsars, as the latter 
have no companion stars and/or disks. But the GPS pulsars appar- 
ently are surrounded by some kind of environment that can affect 
the spectra of those pulsars in the same way as the stellar wind 
affects the B 1259-63 spectrum. As an example we can take the 
cases of B 1823- 13 which is surrounded by a compact pulsar wind 
nebula lPavlov et alj d201 ll), or PSR B 1054-62 wh ich lies behind or 
within a dense HII region dKoribalski et alj 19951) . All GPS pulsars 
have relatively high Dispersion Measures that in some cases, are 
too large to be accounted for by the Galactic electron density, and 
thus, we can speculate that there is a quite high particle density in 
the vicinity of these pulsars (see also lKiiak et alj|201lh . As soon 
as all the necessary observational data are available, we will model 
the physical conditions which cause the evolution of the B 1259-63 
spectrum. Consequently, we will be able to estimate the particle 
density, temperature and magnetic field that are necessary to form 
the spectra shown in Fig. [4] We believe that the same physical pro- 
cesses (i.e. free-free and/or cyclotron absorption) are responsible 
for both B 1259-63 and the GPS pulsars spectra, therefore, we will 
be able to estimate characteristic values for the GPS pulsar sur- 
roundings in a similar way. The only difference could be an in- 
variable shape of the gigahertz-peaked spectra, in contrast with the 
B 1259-63 spectrum, which undergoes evolution due to orbital mo- 
tion. 



5 SUMMARY 

Finally, we can conclude that the B 1259-63 spectrum undergoes 
evolution while the pulsar orbits around the Be star LS 2883. The 
most significant change of the spectrum shape occurs during a cou- 
ple of days just before the eclipse as well as immediately after- 
wards (the thickened black segments in the bottom panel, Fig.[4}. 
We argue that the observed variation of the spectra is caused by a 
combination of two effects: the free-free absorption in the stellar 
wind and the cyclotron resonance in the magnetic field. This field 
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Figure 4. The fits to the B1259-63 spectra for the each orbital phase range (left and right panels). The corresponding data points are given in Fig.[5]and the 
fitted parameters are presented in Table 2. Bottom panel schematically illustrates the orbital epochs over which the spectra are averaged. 



is associated with the disk and is infused by the relativistic parti- 
cles of the pulsar wind. Having noticed the apparent resemblance 
between the B 1259-63 spectrum and the GPS, we suggest that the 
same mechanisms should be responsible for both cases. Thus, the 
case of B 1259-63 can be treated as a key factor to explain the GPS 
phenomenon observed for the solitary pulsars with interesting envi- 
ronments. Therefore, the binary system B1259-63/LS 2883 seems 
to be an important astrophysical laboratory to study the interaction 
between pulsar and t heir environments , such as bow-shocks and 
pulsar wind nebulae JPavlovet all201lh . 
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